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Figure 2.18. Solar radiation on a horizontal surface. 

2.5.3 Solar Radiation on a Tilted Surface 

Solar radiation on an arbitrary tilted surface having a tilt angle of /3 from the horizontal 
and an azimuth angle of a w (assumed + west of south), as shown in Fig. 2.19, is the 
sum of components consisting of beam (I bc ), sky diffuse (I dc ) and ground reflected so¬ 
lar radiation (/ ): 

v r,c' 



(2.42) 


If i is the angle of incidence of the beam radiation on the tilted surface, it is simple to 
show that the instantaneous beam radiation on the surface per unit area is: 


h,c = h.N cos '• 


(2.43) 


From the geometry in Fig. 2.19, it can be shown that the angle of incidence i for the 
surface (angle between the normal to the surface and a line collinear with the sun’s 
rays) is related to the solar angles as 


cos i = cos a cos (a s — tf^sin /3 4- sin a cos /3. 


(2.44) 


The diffuse radiation on the surface (I dc ) can be obtained by multiplying the sky diffuse 
radiation on a horizontal surface by the view factor between the sky and the surface:* 


*The surface has been assumed infinitely large for this view factor. See Section 2.2.3. 
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Incident Angle i 

Figure 2.19. Definitions of solar angles for a tilted surface. 

Uc = 4*0 + cos /3)/2 
= CI bN (\ + cos/3)/2 

= C/ fcJV cos 2 03/2). (2.45) 

The ground reflected solar radiation can be found from the total solar radiation inci¬ 
dent on a horizontal surface and the ground reflectance p as: 

I,c = hP- (2.46) 

The pan of I r intercepted by the tilted surface can be found by multiplying the ground 
reflected radiation by the view factor between the surface* and the ground: 

I rc = p/,,( 1 - cos/3)/2 = pl h sin 2 (/3/2) (2.47) 

= pl b , v (sin a 4- Q sin 2 (/3/2). 

For ordinary ground or grass, p is approximately 0.2 and for snow covered ground it 
can be taken as approximately 0.8. 

Example 2.6(a). Find the instantaneous solar radiation at 12:00 noon Eastern 
Standard Time on a solar collector surface (/3 = 30°. a w = 4-10°) on February 1 in 
Gainesville, FL. 

Solution. From Example 2.1, for February 1: 

n = 32, 8 S = - 17.5°, and ET = - 13.7 min 

*The tilted surface and the ground in front of it have been assumed to be infinitely large for this view 
factor. 
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For finding the values of solar radiation on the collector, we will need to calculate 
angles a , a s , h s , and i. 

Solar Time = LST + ET + 4 (l st — l local ) 

= 12:00 - 13.7 min + 4(75° - 82.27°) 

= 11:17.2 am, 

_ minutes from solar noon 

4 min/deg 

-42.8 

=-= — 10.7° (—before noon). 

4 

From Eq. (2.28), 

a = sin _, (sin(29.68°) sin(- 17.5°) + cos(29.68°) cos(- 17.5°) cos(— 10.7°)) 

= 41.7°. 

From Eq. (2.29), 

a s = sin _1 (cos(- 17.5°)sin(- 10.7°)/cos(41.7°)) 

= -13.7°. 

Angle of incidence i for the solar collector is given by Eq. (2.44): 

cos i = cos(41.7°) cos(— 13.7° - 10°) sin(30°) + sin(41.7°) cos(30°) 
i = 23.4° 

To calculate the solar radiation using Eqs. (2.42, 2.43, 2.45, and 2.47), we need to find 
Kn and L 

Extraterrestrial Solar Radiation 

/ = / 0 [1 + 0.034 cos(360n/365.25)°] 

= 1353W/m 2 [l + 0.034 cos(360x32/365.25)] 

= 1392W/m 2 , 
hjN = C„Ie-V™ a 

(Find k from Table 2.4. Assume C n = 1) 

= 1392e _0144/sin(41 * 7O) 

= 1121 W/m 2 . 
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Beam radiation on the collector (Eq. 2.43): 

4,c = 4 ,n cos 1 

= (1121)cos(23.4°) 

= 1029W/m 2 . 

Sky diffuse radiation on the collector (Eq. 2.45): 

4, = CI bN cos 2 (,8/2) 

= (0.060)(1121)cos 2 (30/2) 

(Find C from Table 2.4.) 


63W/m 2 . 

Ground reflected radiation on the collector (Eq. 2.47): 


4,c = P4^( sin « + C)sin 2 (/3/2) 

Assume 


p = 0.2 

4.c = (0.2)(1121 )(sin41.7° + 0.060)sin 2 (15°) 

= llW/m 2 . 

Total insolation on the collector: 

4 = 1029 + 63 + 11 = 1103W/m 2 . 

Example 2.6(b). Repeat the calculations in Example 2.6(a) for a north facing so¬ 
lar collector (/3 = 30°, a K = 10°) in Buenos Aires, Argentina, (latitude = 34° — 35'S, 
longitude = 58° - 29'W, Standard Meridian = 45°W) 

Solution. Following the procedure in Example 2.6(a), the solar angles and the an¬ 
gle of incidence are found as: 

Solar Time = 12:00 - 13.7 min + 4(45° - 58.48°)min 
= 10:52 am 
h s = -68/4 = -17° 

a = sin-i [sin(—34.58°)sin(—17.5°) + cos(-34.58°)cos(-17.5°)cos(-17 0 )] 

= 67.2° 
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a s = sin- 1 [cos(— 17.5°)sin(— 17°)/cos(67.2°)] 
= -46° 


Therefore, the angle of incidence 

i = cos~ 1 [cos(67.2°)cos(—46° - 10°)sin(30°) + sin(67.2°)cos(30 0 )] 

= 24.95° 

Assuming that the solar radiation model used in Example 2.2(a) applies. (Note: the 
model was developed for locations in the U.S.) 



= 1191 W/m 2 
I bc = 1191 cos(24.95°) 

= 1080 W/m 2 

l dc = (0.060X1191) cos 2 (30°/2) 

= 67 W/m 2 

I rc = (0.2)(1191)(sin 67.2° + 0.060) sin 2 (15°) 
= 16 W/m 2 


Total insolation on the collector 


/ = 1080 + 67 +16 


= 1163 W/m 2 


2.5.4 Monthly Solar Radiation Estimation Models 

One of the earliest methods of estimating solar radiation on a horizontal surface was 
proposed by the pioneer spectroscopist Angstrom. It was a simple linear model relating 
average horizontal radiation to clear-day radiation and to the sunshine level, that is, 
percent of possible hours of sunshine. Since the definition of a clear day is somewhat 
nebulous, Page [57] refined the method and based it on extraterrestrial radiation in¬ 
stead of the ill-defined clear day: 




(2.48) 
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where H h and H oh are the horizontal terrestrial, and horizontal extraterrestrial radiation 
levels averaged for a month, PS is the monthly averaged percent of possible sunshine 
(that is, hours of sunshine/maximum possible duration of sunshine X 100), a and b are 
constants for a given site, and n and N are the monthly average numbers of hours of 
bright sunshine and day length respectively. The ratio n/N is also equivalent to the 
monthly average percent sunshine {PS). H o h can be calculated by finding H o h from the 
following equation (2.49), using Eqs. (2.35) and (2.44) and averaging I o h for the num¬ 
ber of days in each month; or data in Appendix 2 (Table A2.2) can be used. 



(2.49) 


Some typical values of a and b are given in Table 2.5 [46]. Additional values for 
worldwide locations are given in Appendix 2 (Table A2.4). 

Example 2.7. Using the predictive method of Angstrom or Page, estimate the 
monthly solar radiation for the North Central Sahara Desert (Tamanrasset, Algeria 
area) at latitude = 25°N. Percentages of possible sunshine and extraterrestrial radia¬ 
tion for this site are given in the table below. 


Solution. Using the climate data given, the expected monthly average horizontal 
radiation for the North Sahara is calculated in the following table using a = 0.30 and 
b = 0.43 from Table 2.5. 










Ho* 


«k 


Month 

PS/100 

(kJ/m 2 • day) 

kJ/m 2 • day 


Btu/ft 2 • day 

Jan. 

0.88 

23,902 

16,215 


1425 

Feb. 

0.83 

28,115 

18,469 


1626 

Mar. 

0.90 

32,848 

22,143 


1950 

Apr. 

0.85 

37,111 

24,697 


2174 

May 

0.80 

39.356 

25,345 


2231 

June 

0.76 

40,046 

25,101 


2210 

Jul. 

0.86 

39,606 

36,528 


2336 

Aug. 

0.83 

37,832 

24,852 


2188 

Sep. 

0.77 

34,238 

21,608 


1902 

Oct. 

0.86 

29,413 

19,701 


1735 

Nov. 

0.85 

24,909 

16,577 


1460 

Dec. 

0.80 

22,669 

14,599 


1285 


“Monthly averaged, daily extraterrestrial radiation. 


The high levels of radiation predicted above are typical of the North Sahel region and 
are higher than most U.S. locations except the Mojave Desert. 

A number of researchers found Angstrom-Page type correlations for specific loca¬ 
tions which are listed in Table 2.6. Some of these include additional parameters such as 
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Table 2.5. Coefficients a and b in the Angstrom-Page regression equation* 

Sunshine hours in 
percentage of possible 


Location 

Climate^ 

Range 

Avg. 

a 

b 

Albuquerque. NM 

BS-BW 

68-85 

78 

0.41 

0.37 

Atlanta, GA 

Cf 

45-71 

59 

0.38 

0.26 

Blue Hill. MA 

Df 

42-60 

52 

0.22 

0.50 

Brownsville, TX 

BS 

47-80 

62 

0.35 

0.31 

Buenos Aires, Arg. 

Cf 

47-68 

59 

0.26 

0.50 

Charleston, SC 

Cf 

60-75 

67 

0.48 

0.09 

Dairen, Manchuria 

Dw 

55-81 

67 

0.36 

0.23 

El Paso, TX 

BW 

78-88 

84 

0.54 

0.20 

Ely, NV 

BW 

61-89 

77 

0.54 

0.18 

Hamburg, Germany 

Cf 

11-49 

36 

0.22 

0.57 

Honolulu, HI 

Af 

57-77 

65 

0.14 

0.73 

Madison, WI 

Df 

40-72 

58 

0.30 

0.34 

Malange. Angola 

Aw-BS 

41-84 

58 

0.34 

0.34 

Miami, FL 

Aw 

56-71 

65 

0.42 

0.22 

Nice, France 

Cs 

49-76 

61 

0.17 

0.63 

Poona. India (monsoon) 

Am 

25—49 

37 

0.30 

0.51 

(dry) 


65-89 

81 

0.41 

0.34 

Stanleyville, Congo 

Af 

34-56 

48 

0.28 

0.39 

Tamanrasset, Algeria 

BW 

76-88 

83 

0.30 

0.43 


+ From Lof et al. [46] with permission. 

++ Af = tropical forest climate, constantly moist, rainfall all through the year 

Am = tropical forest climate, monsoon rain, short dry season, but total rainfall sufficient to support 
rain forest 

Aw = tropical forest climate, dry season in winter 
BS = steppe or semiarid climate 
BW = desert or arid climate 

Cf = mesothermal forest climate, constantly moist, rainfall all through the year 
Cs = mesothermal forest climate, dry season in winter 

Df = microthermal snow forest climate, constantly moist, rainfall all through the year 
Dw = microthermal snow forest climate, dry season in winter 


relative humidity and ambient temperature. Correlations listed in the table may be 
used for the specific locations for which they were developed. 

Another meteorological variable that could be used for solar radiation prediction 
is the opaque cloud cover recorded at many weather stations around the world. This 
quantity is a measure of the percent of the sky dome obscured by opaque clouds. Be¬ 
cause this parameter contains even less solar information than sunshine values, it has 
not been useful in predicting long-term solar radiation values. A subsequent section, 
however, will show that cloud cover, when used with solar altitude angle or air mass, is 
a useful estimator of hourly direct radiation. 
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Table 2.6. Angstrom-Page type correlations for specific locations 



Measured data 



Authors 

correlated 

Correlation equations 

* 

Iqbal [28] 

Canada, 

3 locations 

= 0.791 - 0.635 

H i, 

( 1 ) 



S* = 0.163 + 0.4781 

H h ' 




-0.176 + ,.«(|)-... 2 (|)’ 


Garg [17] 


India, 

11 locations, 
20 years’ data 




H 


ojl 


= 0.3156 + 0.45201 — 


Dl 


H, 


= 0.3616 - 0.2123 


ojl 


° h = 0.8677 - 0.7365( = 


H, 


Hussain [25] India 


Coppolino [11] Italy 


3_ 

Hoj, 

2k. 

H /«\ 045 

sin( ““ )o ° 5 

a sn = Solar Elevation at noon on the 15 th of each month, 
degrees 

0.15 =£ = £ 0.90 
N 


= 0.394 + 0.3641 ^ 7 1 - 0.0035 W„, 
= 0.306 — 0.165 =, — 0.0025 W,,, 


Akinoglu & Italy 

Ecevit [2] 


^-0 |45 *0.»5(|)-0.I80(|) : 


©■ 02M+o ' 758 (I)-°- 25o {[(I)T + 1) 

0.326(|) - 0.433(f)’ 


an = 0.035 + 

N 


Ogelman et al. [56] 


Turkey, 

2 locations, 

3 years’ data 
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Table 2.6. ( Continued) 


Authors 


Measured data 
correlated 


Correlation equations* 


Gopinathan [ 19] 


40 locations 



= a + b 



around the 
world 



a = - .309 + .539 cos L- .0639 h + 0.290 


b = 1.527 - 1.027 cos L + .0926 h -.359 


*H a , H h , H o h , D h are monthly averaged daily values. 

N' = maximum duration for which Campbell-Stokes recorder can be active, i.e.. Solar Elevation >5°. 
W (it = relative humidity X (4.7923 + 0.3647T 4- 0.0557^ -I- 0.00037 3 
T = ambient temperature °C 
W at = gm moisture/m 3 
h = Elevation in km above sea level. 

L = latitude 

2.6 MODELS BASED ON LONG-TERM MEASURED HORIZONTAL 
SOLAR RADIATION 

Long-term measured solar radiation data is usually available as monthly averaged total 
solar radiation per day on horizontal surfaces. In order to use this data for tilted sur¬ 
faces, the total solar radiation on a horizontal surface must first be broken down into 
beam and diffuse components. A number of researchers have proposed models to do 
that, prominent among them being Liu and Jordan, Collares-Pereira and Rabl, and 
Erbs, Duffie and Klein. 

2.6.1 Monthly Solar Radiation on Tilted Surfaces 

In a series of papers, Liu and Jordan [41-45] have developed an essential simplifica¬ 
tion in the basically complex computational method required to calculate long-term ra¬ 
diation on tilted surfaces. This is called the LJ method. The fundamental problem in 
such calculations is the decomposition of long-term measured total horizontal radia¬ 
tion into its beam and diffuse components. 

If the decomposition can be computed, the trigonometric analysis presented ear¬ 
lier can be used to calculate incident radiation on any surface in a straightforward man¬ 
ner. Liu and Jordan (LJ) correlated the diffuse-to-total radiation ratio ( DjH h ) with the 
monthly clearness index K T , which is defined as 


Ok 


(2.50) 



where H h is the monthly averaged terrestrial radiation per day on a horizontal surface. 
H oh is the corresponding extraterrestrial radiation, which can be calculated from Eq. 
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(2.49) by averaging each daily total for a month. The original LJ method was based 
upon the extraterrestrial radiation at midmonth, which is not truly an average. 

The LJ correlation predicts the monthly diffuse ( D h ) to monthly total H h ratio. It 
can be expressed by the empirical Eq. 

= 1.390 - 4.027 K t + 5.531^| - 3.108^ . (2.51) 

H h 


Note that the LJ correlation is based upon a solar constant value of 1394 W/m 2 
(442 Btu/hr • ft 2 ), which was obtained from terrestrial observations, whereas the 
newer value, based on satellite data, is 1377 W/m 2 (437 Btu/hr • ft 2 ). The values of 
K T must be based on this earlier value of the solar constant to use the LJ method. Col- 
lares-Pereira and Rabl [10] conducted a study and concluded that although Liu and Jor¬ 
dan’s approach is valid, their correlations would predict significantly smaller diffuse ra¬ 
diation components. They also concluded that Liu and Jordan were able to correlate 
their model with the measured data because they used the measured data which was not 
corrected for the shade ring (see solar radiation measurements). Collares-Pereira and 
Rabl (C-P&R) also introduced the sunset hour angle h ss in their correlation to account 
for the seasonal variation in the diffuse component. The C-P&R correlation is: 


= 0.775 + 0.347^ 



0.505 + 0.0261 




cos(2AT r - 1.8), 

(2.52) 


where h ss is the sunset hour angle in radians. The C-P&R correlation agrees well with 
the correlations for India [8], Israel [71] and Canada [66] and is, therefore, preferred to 
equation (2.51). 

The monthly average beam component B h on a horizontal surface can be readily 
calculated by simple subtraction since D h is known: 

B h = H h - D h . (2.53) 

It will be recalled on an instantaneous basis from Eqs. (2.41) and (2.43) and Fig. (2.19) 
that 


b,h 


1 h,N sin a 

4,C = 4.v cos *•» 


(2.54) 

(2.43) 


where l bh is the instantaneous horizontal beam radiation. Solving for I bc , the beam ra¬ 
diation on a surface, 


cos l 


sin a 


4,c — 4./i 


(2.55) 






FUNDAMENTALS OF SOLAR RADIATION 


57 


The ratio in parentheses is usually called the beam radiation tilt factor R b . It is a purely 
geometric quantity that converts instantaneous horizontal beam radiation to beam radi¬ 
ation intercepted by a tilted surface. 

Equation (2.55) cannot be used directly for the long-term beam radiation B h . To be 
strictly correct, the instantaneous tilt factor R b should be integrated over a month with 
the beam component I hh used as a weighting factor to calculate the beam tilt factor. 
However, the LJ method is used precisely when such short-term data as I bh are not 
available. The LJ recommendation for the monthly mean tilt factor R b is simply to cal¬ 
culate the monthly average of cos i and divide it by the same average of sin a. In equa¬ 
tion form for south-facing surfaces, this operation yields: 


Ru = 


cos(L — /3)cos S s sin h sr 4- h sr sin(L - /3) sin 8 S 
cos L cosS v sin h sr (a = 0) + hJa = 0) sin L sin 5 ’ 


(2.56) 


where the sunrise hour angle h gr (a = 0) in radians is given by Eq (2.30) and h sr is the 
min [I h s (a = 0)1, I h s (i = 90°)l], respectively, and are evaluated at midmonth. Non¬ 
south-facing surfaces require numerical integration or iterative methods to determine 
R b . The long-term beam radiation on a tilted surface B c is then, 




(2.57) 


which is the long-term analog of Eq. (2.43). Values of R h are tabulated in Appendix 2, 
Table A2.5. 

Diffuse radiation intercepted by a tilted surface differs from that on a horizontal 
surface, because a tilted surface does not view the entire sky dome, which is the source 
of diffuse radiation. If the sky is assumed to be an isotropic source of diffuse radiation, 
the instantaneous and long-term tilt factors for diffuse radiation, R d and R d respec¬ 
tively, are equal and are simply the radiation view factor from the plane to the visible 
portion of a hemisphere. In equation form: 


R d = R d = cos 2 — = (1 + cos/3)/2. 


(2.58) 


In some cases where solar collectors are mounted near the ground, some beam and 
diffuse radiation reflected from the ground can be intercepted by the collector surface. 
The tilt factor R r for reflected total radiation (D h + B h ) is then calculated to be 


R m = 


R 




n , 0 = psin^t- = p(l - cos/3)/2, 

Ui. -r Du l 


(2.59) 


in which p is the diffuse reflectance of the surface south of the collector assumed uni¬ 
form and of infinite extent. 

For snow, p = 0.75; for grass and concrete, p = 0.2. A more complete list of re¬ 
flectances is provided in Table A2.7 of Appendix 2. The total long-term radiation inter- 
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cepted by a surface H c is then the total of beam, diffuse, and diffusely reflected compo¬ 
nents: 


(2.60) 


H e = R b B k + R 4 D k + RJtP k + B k ). 


Using Eqs. (2.61) and (2.62), we have 



(2.61) 


in which R b is calculated from Eq. (2.56). 

Example 2.8. Using the H h data calculated in Example 2.7 in place of the long¬ 
term measured data for the North Central Sahara Desert at latitude 25°N, find the 
monthly averaged insolation per day on a south-facing solar collector tilted at an angle 
of 25° from the horizontal. 

Solution. The solution below is for the month of January. Values for the other 
months can be found by following the same method. 


H h = 16,215kJ/m 2 - day. 


From Table A2.2a: 


H oJl = 23,902. 


Therefore, 


= HjHo* = 0.678. 


8 s and h sr can be found for the middle of the month (January 16). 


8 S = 23.45° sin[360(284 + 16)365]° 

= -21.1°, 


h sr {a = 0) = -cos ‘(-tan L tan 8) 


= —79.6° or — 1.389rad. 


and/? = 1.389 

55 


Using CP & R correlation. 
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= 0.775 + 0.347(1.389 - 1.5708) 

H h 

- [0.505 + 0.0261(1.389 - 1,5708)]cos(2 X 0.678 - 1.8) 

= 0.212 

Therefore, 

D h = 0.212 X 16,215 = 3,438kJ/m 2 - day 
and B h = H h - D h = 12,777kJ/m 2 - day 

Insolation on a tilted surface can be found from Eq. (2.60). We need to find R B from 
Eq. (2.56). 

Therefore, 

- _ cos(O) cos(—21.1°) sin( —79.6°) - 1.389sin(0) sin(-21.1°) 

Rb ~ cos(25°) cos(—21.1°) sin(—79.6°) - 1.389sin(25°) sin(-21.1°) 

= 147. 

R d = cos 2 (25/2) = 0.953. 

R r = p sin 2 (/3/2) (Assume p = 0.2) 

= 0.2 sin 2 (12.5°) 

= 0.009. 

Therefore, 

H c = (1.47)(12,777) + 0.953(3,438) + 0.009(16,215) 

= 22,205 kJ/m 2 . 


2.6.2 Circumsolar or Anisotropic Diffuse Solar Radiation 

The models described in the above sections assume that the sky diffuse radiation is 
isotropic. However, this assumption is not true because of circumsolar radiation 
(brightening around the solar disk). Although the assumption of isotropic diffuse solar 
radiation does not introduce errors in the diffuse values on horizontal surfaces, it can 
result in errors of 10 to 40% in the diffuse values on tilted surfaces. A number of re¬ 
searchers have studied the anisotropy of the diffuse solar radiation because of circum¬ 
solar radiation. Temps and Coulson [73] introduced an anisotropic diffuse radiation al¬ 
gorithm for tilted surfaces for clear sky conditions. Klucher [37] refined the Temps and 
Coulson algorithm by adding a cloudiness function to it: 
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R d = —(1 + cos /3)M,M 2 , 


(2.62) 


where 


M ] = 1 + F sin 3 (/3/2), 

M 2 = 1 + F cos 2 i sin 3 (z), and 

f= i - (ZVtf*) 2 - 


(2.63) 

(2.64) 

(2.65) 


Examining F, we find that under overcast skies ( D h = H h ), R (J in Eq. 2.62 reduces 
to the isotropic term of Liu and Jordan. The Klutcher algorithm reduces the error in 
diffuse radiation to about 5%. 

In summary, monthly averaged, daily solar radiation on a surface is calculated by 
first decomposing total horizontal radiation into its beam and diffuse components us¬ 
ing Eq. (2.51) or (2.52). Various tilt factors are then used to convert these horizontal 
components to components on the surface of interest. 


2.6.3 Daily Solar Radiation on Tilted Surfaces 

Prediction of daily horizontal total solar radiation for sites where solar data are not 
measured can be done using the Angstrom-Page model. Instead of monthly values, 
however, daily values are used for percent sunshine PS and extraterrestrial radiation 
I day . The results of using this simple model would be expected to show more scatter 
than monthly values, however. 

All U.S. National Weather Service (NWS) stations with solar capability report 
daily horizontal total (beam and diffuse) radiaton. Liu and Jordan have extended their 
monthly method described above to apply to daily data. The equation, analogous to 
Eq. (2.51), used to calculate the daily diffuse component I dh is [6] 


kb. = i.oo 45 + 0.04349AV - 3.5227 K T 2 + 2.6313 K T \ and (2.66) 


h 


K T < 0.75, 


where K T (no overbar) is the daily clearness index analogous to the monthly K r . (In 
this section, overbars indicate daily radiation totals.) For values of K T > 0.75, the dif- 
fuse-to-total ratio is constant at a value of 0.166. K T is given by 



(2.67) 


The daily extraterrestrial total radiation I oh is calculated from Eq. (2.49). Note that Eq. 
(2.66) is based on the early solar constant value of 1394 W/m 2 (442 Btu/hr • ft 2 ). 
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The daily horizontal beam component I hh is given by simple subtraction: 

= h~ W ( 2 . 68 ) 

The beam, diffuse, and reflected components of radiation can each be multiplied by 
their tilt factors R ., R „ and R to calculate the total radiation on a tilted surface 

n a 7 r 

l = RJ bJ , + RJ dM + R,0b.h + /*»). ( 2 - 69 > 


in which 


cos(L — /3)cos3 5 sin h sr + /? 5r sin(L — /3)sin^_ 
cos L cos6 5 sin h sr (a = 0) + h sr (a = 0)sin L sin5/ 





, and 


(2.70) 

(2.71) 


/? 


r 


P 



(2.72) 


by analogy with the previous, monthly analysis. 

If daily solar data are available, they can be used for design, the same as monthly 
data. Daily calculations are necessary when finer time-scale performance is required. 
In addition, daily data can be decomposed into hourly data, which are useful for calcu¬ 
lations made with large, computerized solar system simulation models. 


2.6.4 Hourly Solar Radiation on Tilted Surfaces 

Hourly solar radiation can be predicted in several ways. Correlations between hourly 
total and hourly diffuse (or beam) radiation or meteorological parameters such as 
cloud cover or air mass may be used. Alternatively, a method proposed by Liu and Jor¬ 
dan based on the disaggregation of daily data into hourly data could be used [41-45]. 
Even if hourly NWS data are available, it is necessary to decompose these total values 
into beam and diffuse components depending upon the response of the solar conver¬ 
sion device to be used to these two fundamentally different radiation types. 

Randall and Leonard [62] have correlated historical data from the NWS stations at 
Blue Hill, Massachusetts, and Albuquerque, New Mexico to predict hourly beam radi¬ 
ation I b and hourly total horizontal radiation I h . This method can be used to decom¬ 
pose NWS data into its beam and diffuse components. 

The hourly beam radiation was found to be fairly well correlated by hourly per¬ 
cent of possible insolation k r defined as 



(2.73) 
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in which I gh is the hourly horizontal extraterrestrial radiation, which can be evaluated 
from Eq. (2.49) using one-hour integration periods. Carrying out the integration over a 
one hour period yields 


K,h ~ I + 0.034 cos j(0.9972 cos L coscos h s -I- sin L sin S^), (2.74) 

where the solar hour angle h s is evaluated at the center of the hour of interest. 

The direct normal-beam correlation based on k t is [45] 

hjt = -520 + 1800 k,(W/m 2 ) 0.85 > k, > 0.30 (2.75) 

h.N = 0 k,< 0.30 (2.76) 

This fairly simple correlation gives more accurate I bN values than the more cumber¬ 
some Liu and Jordan procedure, at least for Blue Hill and Albuquerque. Vant-Hull and 
Easton [80] have also devised an accurate predictive method for beam radiation. 

Randall and Leonard [62] have made a correlation of total horizontal hourly radi¬ 
ation l b (W/m 2 ) on the basis of opaque cloud cover CC and air mass m using data for 
Riverside, Los Angeles, and Santa Monica, California. Cloud cover is defined as CC = 

1.0 for fully overcast and CC = 0.0 for clear skies. A polynomial fit was used: 

l nh 

K = TOO (83 '° 2 ~ 3M7m ~ 4.407CC + 1.1013CC 2 - 0.1109CC 3 ) (2.77) 

The average predictive error for this correlation was ± 2.3 percent of the NWS data; 
the correlation coefficient of I h with CC is 0.76 for the data used. Equation (2.80) was 
used to predict I h for Inyokern, California, a site not used in the original correlation. 
Predictions of solar radiation for Inyokern were within 3.2 percent of NWS. The dif¬ 
fuse radiation can be calculated from Eqs. (2.75), (2.76), and (2.77): 

hb = h ~ 4sina (2.78) 

The American Society of Heating, Refrigerating, and Air Conditioning Engineers 
(ASHRAE) has calculated hourly clear-sky values on vertical, horizontal, and tilted 
surfaces. One table for a full year has been prepared for each of six values of 
latitude spanning the continental United States. In addition, tables of solar azimuth and 
altitude angle have been prepared for the same six latitudes. The solar radiation 
data represent the maximum that could be expected on a clear day and are therefore of 
limited usefulness in the design of solar systems. They can be used, however, to cal¬ 
culate cooling loads on buildings and the like. The ASHRAE tables are contained in 
Appendix 2. 
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supersaturated fluosilicic acid solution, the acid attacks the glass and leaves a porous 
silica surface layer. This layer has an index of refraction intermediate between glass and 
air. Figure 3.5 shows the spectral reflectance of a pane of glass before and after etching. 


3.2 FLAT PLATE COLLECTORS 

A simple flat plate collector consists of an absorber surface (usually a dark, thermally 
conducting surface); a trap for reradiation losses from the absorber surface (such as 
glass which transmits shorter wavelength solar radiation but blocks the longer wave¬ 
length radiation from the absorber); a heat transfer medium such as air, water, etc.; and 
some thermal insulation behind the absorber surface. Flat plate collectors are used typ¬ 
ically for temperature requirements up to 75°C although higher temperatures can be 
obtained from high efficiency collectors. These collectors are of two basic types based 
on the heat transfer fluid: 

Liquid type: where heat transfer fluid may be water, mixture of water and antifreeze, 
oil, etc.; 

Air type: where heat transfer medium is air (used mainly for drying and space heating 
requirements). 

3.2.1 Liquid-Type Collectors 

Figure 3.6 shows a typical liquid-type flat plate collector. In general it consists of: 

1. Glazing : one or more covers of transparent material like glass, plastics, etc. Glaz¬ 
ing may be left out for some low-temperature applications. 



ABSORBER PLATE 


INSULATION 


Figure 3.6. Typical liquid-type flat plate collector (courtesy of Morning Star Corporation, 
Orange Park. Florida). 


GLAZING 
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2. Absorber : a plate with tubes or passages attached to it for the passage of a working 
fluid. The absorber plate is usually painted flat black or electroplated with a selec¬ 
tive absorber. 

3. Headers or manifolds: to facilitate the flow of heat transfer fluid. 

4. Insulation : to minimize heat loss from the back and the sides. 

5. Container, box or casing. 


3.2.2 Air-Type Collectors 

Air types of collectors are more commonly used for agricultural drying and space heat¬ 
ing applications. Their basic advantages are low sensitivity to leakage and no need for 
an additional heat exchanger for drying and space heating applications. However, be¬ 
cause of the low heat capacity of the air and the low convection heat transfer coeffi¬ 
cient between the absorber and the air, a larger heat transfer area and higher flow rates 
are needed. Figure 3.7 shows some common configurations of air heating collectors. 
Common absorber materials include corrugated aluminum or galvanized steel sheets, 
black metallic screens, or simply any black painted surface. 

Unglazed, transpired solar air collectors offer a low-cost opportunity for some ap¬ 
plications such as preheating of ventilation air and agricultural drying and curing [49]. 
Such collectors consist of perforated absorber sheets that are exposed to the sun and 
through which air is drawn. The perforated absorber sheets are attached to the vertical 
walls, which are exposed to the sun. Kutcher and Christensen [50] have given a de- 
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Figure 3.7. Some common configurations of air heating collectors. 
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tailed thermal analysis of unglazed transpired solar collectors. (See Chapter 5, Section 
5.4.5 for additional details.) 

The most important components, whose properties determine the efficiency of so¬ 
lar thermal collectors, are glazings and absorbers. 

3.2.3 Glazings 

The purpose of a glazing or transparent cover is to transmit the shorter wavelength so¬ 
lar radiation but block the longer wavelength reradiation from the absorber plate, and 
to reduce the heat loss by convection from the top of the absorber plate. Consequently, 
an understanding of the process and laws that govern the transmission of radiation 
through a transparent medium is important. Section 3.1 describes in brief the transmis¬ 
sion of radiation through materials. 

Glass is the most common glazing material. Figure 3.8 shows transmittance of glass 
as a function of wavelength. Transparent plastics, such as polycarbonates and acrylics 
are also used as glazings for flat plate collectors. The main disadvantage of plastics is 
that their transmittance in the longer wavelength is also high, therefore, they are not as 
good a trap as glass. Other disadvantages include deterioration over a period of time due 
to ultraviolet solar radiation. Their main advantage is resistance to breakage. Although 
glass can break easily, this disadvantage can be minimized by using tempered glass. 

In order to minimize the upward heat loss from the collector, more than one trans¬ 
parent glazing may be used. However, with the increase in the number of cover plates, 
transmittance is decreased. Figure 3.9 shows the effect of number of glass cover plates 
on transmittance. 



VW 






92 


PRINCIPLES OF SOLAR ENGINEERING 



Angle of Incidence 

Figure 3.9. Transmittance of multiple glass covers having an 
index of refraction 1.526". 


Absorbers. The purpose of the absorber is to absorb as much of the incident solar ra¬ 
diation as possible, re-emit as little as possible, and allow efficient transfer of heat to a 
working fluid. The most common forms of absorber plates in use are shown in Figure 
3.10. The materials used for absorber plates include copper, aluminum, stainless steel, 
galvanized steel, plastics and rubbers. Copper seems to be the most common material 
used for absorber plates and tubes because of its high thermal conductivity and high 
corrosion resistance. However, copper is quite expensive. For low-temperature appli¬ 
cations (up to about 50°C or 120°F) a plastic material called ethylene propylene poly¬ 
mer (trade names EPDM, HCP, etc.) can be used to provide inexpensive absorber ma- 





Roll-Bond 

Figure 3.10. Common types of absorber plates. 
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HCP® Absorber Panel 


Manifold 


Frame 


Substrate 


Back Brace 


Figure 3.11. Cutaway view of a typical collector made from ethylene propylene polymer (Courtesy of 
Sealed Air Corporation, Saddlebrook, NJ). 


terial. To compensate for the low thermal conductivity, a large surface area is provided 
for heat transfer. Figure 3.11 shows a typical collector made from such material. 

In order to increase the absorption of solar radiation and to reduce the emission 
from the absorber, the metallic absorber surfaces are painted or coated with flat black 
paint or some selective coating. A selective coating has high absorptivity in the solar 
wavelength range (0.3 to 3.0 ^m). Absorptivities and emissivities of some common se¬ 
lective surfaces are given in Table 3.3. 

A simple and inexpensive collector consists of a black painted corrugated metal 
absorber on which water flows down open, rather than enclosed in tubes. This type of 
collector is called a trickle collector and is usually built on-site. Although such a col¬ 
lector is simple and inexpensive, it has the disadvantages of condensation on the glaz¬ 
ing and a higher pumping power requirement. 

3.2.4 Energy balance for a flat-plate collector. The thermal performance of any type 
of solar thermal collector can be evaluated by an energy balance that determines the 
portion of the incoming radiation delivered as useful energy to the working fluid. For a 
flat-plate collector of an area A c this energy balance on the absorber plate is 

= Q u + <hoss + < 3 - 8 ) 


where 1 = solar irradiation on a collector surface, 

c 

r = effective solar transmittance of the collector cover(s), 

s 

a s = solar absorptance of the collector-absorber plate surface, 
q u = rate of heat transfer from the collector-absorber plate to the working fluid, 
q = rate of heat transfer (or heat loss) from the collector-absorber plate to the 
surroundings, and 

de c /dt = rate of internal energy storage in the collector. 

The instantaneous efficiency of a collector rj c is simply the ratio of the useful energy 
delivered to the total incoming solar energy, or 
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(3.9) 


In practice, the efficiency must be measured over a finite time period. In a standard 
performance test, this period is on the order of 15 or 20 min, whereas for design, the 
performance over a day or over some longer period t is important. Then we have for 
the average efficiency 



(3.10) 


where t is the time period over which the performance is averaged. 

A detailed and precise analysis of the efficiency of a solar collector is complicated 
by the nonlinear behavior of radiation heat transfer. However, a simple linearized 
analysis is usually sufficiently accurate in practice. In addition, the simplified analyti¬ 
cal procedure is very important because it illustrates the parameters of significance for 
a solar collector and how these parameters interact. For a proper analysis and interpre¬ 
tation of these test results an understanding of the thermal analysis is imperative, al¬ 
though for design and economic evaluation the results of standardized performance 
tests are generally used. 


Working fluid 



Tube 


Figure 3.12. Schematic diagram of solar collector. 
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Collector heat-loss conductance. In order to obtain an understanding of the parame¬ 
ters determining the thermal-efficiency of a solar collector, it is important to develop 
the concept of collector heat-loss conductance. Once the collector heat-loss conduc¬ 
tance U c is known, and when the collector plate is at an average temperature 7\, the 
collector heat loss can be written in the simple form 

qioss=U c A c (T c -T a ) (3.11) 

The simplicity of this relation is somewhat misleading because the collector heat-loss 
conductance cannot be specified without a detailed analysis of all the heat losses. Fig¬ 
ure 3.12 shows a schematic diagram of a single-glazed collector, while Fig. 3.13(a) 
shows the thermal circuit with all the elements that must be analyzed before they can 
be combined into a single conductance element shown in Figure 3.13(b). The analysis 
below shows an example of how this combination is accomplished. 

In order to construct a model suitable for a thermal analysis of a flat-plate collec¬ 
tor, the following simplifying assumptions will be made: 




Figure 3.13. Thermal circuits for a flat-plate collector shown in Fig. 3.12: (a) detailed circuit; (b) 
approximate, equivalent circuit to (a). In both circuits, the absorber energy is equal to a s I s , where I s = rl c . 
Collector assumed to be at uniform temperature 7\ 












96 


PRINCIPLES OF SOLAR ENGINEERING 


1. The collector is thermally in steady state. 

2. The temperature drop between the top and bottom of the absorber plate is negligible. 

3. Heat flow is one-dimensional through the cover as well as through the back in¬ 
sulation. 

4. The headers connecting the tubes cover only a small area of the collector and pro¬ 
vide uniform flow to the tubes. 

5. The sky can be treated as though it were a black-body source for infrared radiation 
at an equivalent sky temperature. 

6 . The irradiation on the collector plate is uniform. 

For a quantitative analysis let the plate temperature be T c and assume solar energy 
is absorbed at the rate / a y Part of this energy is then transferred as heat to the working 
fluid, and if the collector is in the steady state, the other part is lost as heat to the ambi¬ 
ent air if T c > 7\ Some of the heat loss occurs through the bottom of the collector. It 
passes first through the back to the environment. Since the collector is in steady state, 
according to Eq. (3.8), 


Qu floss'* (3.12) 

where q hss can be determined using the equivalent thermal circuit as shown in figure 3.13. 

q loss = U C A C {T C - TJ = T “\ (3.13) 

There are three parallel paths to heat loss from the hot collector absorber plate at T to 
the ambient at T a : the top, bottom and edges. Because the edge losses are quite small 
compared to the top and the bottom losses, they are quite often neglected. However, 
they can be estimated easily if the insulation around the edges is of the same thickness 
as the back. The edge loss can be accounted for by simply adding the areas of the back 
(A c ) and the edges (A e ) for back heat loss. 

Therefore, the overall heat loss coefficient is: 


U C A C 



The thermal resistances can be found 


A c + A c + A 

Ry + R 2 , R lb + /? 2 * 
easily from the definition. For example. 


(3.14) 


= ~T and R 2b = 7 — , (3.15) 

*1 "'c,bottom 

where k. and l. are, respectively, the thermal conductivity and thickness of the insula¬ 
tion, and h cbottom is the convective heat transfer coefficient between the collector and 
the air below the collector. In a well-insulated collector, R 2b is much smaller than R U) 
and usually neglected. Referring to figure 3.12 
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A e = 2(/, + / 2 )/ 3 (3.16) 

Since the heat loss from the top is by convection and radiation, it is more compli¬ 
cated than the bottom heat loss. Convection and radiation provide two parallel paths 
for heat loss from the absorber plate at T c to the glass cover at T g , and from the glass 
cover to the ambient. That is, the series resistance of R u and R 2t consists of: 



Therefore, 


1 1 1 

— = — 4- — - h r i + h r , and 

p p p r » 1 c * 1 

'V.l K c.\ 


ill,, 

— =- 1 -= h + h. 

R„ R. 


r, w c,t» 


(3.17) 

(3.18) 


Since thermal radiative heat transfer is proportional to the fourth power of the 
temperature, R r and h r are found as follows: 

Radiative heat transfer from the plate to the glass cover 


4r 


= crA, 


(T 4 C 


T P 


( 1 / e P .i + Ue gJ - 1 ) 


= h rA A c (T c - T g ), 


(3.19) 


where 


6 . = infrared emittance of the plate 

P'i r 

6 . = infrared emittance of the glass cover. 

g’t 

Therefore, 


o<r c + TjiJl + rp 

(i/v + 1 / € g .i ~ i) 


(3.20) 


Similarly, from the radiative heat transfer between the glass plate (at T ? ) and the sky (at 
T sk ) we can find that: 


~ T sky) = h r.~ A c( T g ~ T a)’ 


( 3 . 21 ) 
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or 


^r,oo ^tky)/(T g T a ). (3.22) 

Evaluation of the collector heat-loss conductance defined by Eq. (3.14) requires 
iterative solution of Eqs. (3.19) and (3.21), because the unit radiation conductances are 
functions of the cover and plate temperatures, which are not known a priori. A simpli¬ 
fied procedure for calculating U c for collectors with all covers of the same material, 
which is often sufficiently accurate and more convenient to use, has been suggested by 
Hottel and Woertz [36] and Klein [41]. It is also suitable for application to collectors 
with selective surfaces. For this approach the collector top loss in watts is written in 
the form [1]: 


* = _ Vc ~ T Mr _ 

^ n /(C/T c )[(T c - T a )/(N + f)] 033 + 1 /h^ 

+ ___ o{T*-Tl)A c 

\/[e pJ + 0.05M1 - £ pi< )] + (2N+f- \)/ £gJ - N ' (3 ' 23) 

where/ = (1 - 0.04 h cx + 0.005/i c J)(l + 0.091 N), 

C = 250[1 - 0.004403 - 90)j, 

N = number of covers, 
h cx = 5.7 + 3.8V, and 
6 s .j = infrared emittance of the covers, 

V = wind speed in m/sec. 

The values of q, oploss calculated from Eq. (3.23) agreed closely with the values ob¬ 
tained from Eq. (3.22) for 972 different observations encompassing the following con¬ 
ditions: 


320 <T. <420 K 
260 <T a <310 K 
0.1 <ej <0.95 
0sV< 10 m/sec 
1SN<3 
0 < p < 90 


The standard deviation of the differences in U c = q, low /A c (T t -T a ) was 0.14 W/m^ • 
K for these comparisons. 


3.2.5 Thermal Analysis of Flat-Plate Collector-Absorber Plate 

In order to determine the efficiency of a solar collector, the rate of heat transfer to the 
working fluid must be calculated. If transient effects are neglected [35,41,79], the rate 
of heat transfer to the fluid flowing through a collector depends on the temperature of 
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the collector surface from which heat is transferred by convection to the fluid, the tem¬ 
perature of the fluid, and the heat-transfer coefficient between the collector and the 
fluid. To analyze the rate of heat transfer consider first the condition at a cross section 
of the collector with flow ducts of rectangular cross sections as shown in Fig. 3.14. So¬ 
lar radiant energy impinges on the upper face of the collector plate. A part of the total 
solar radiation falls on the upper surface of the flow channels, while another part is in¬ 
cident on the plates connecting any two adjacent flow channels. The latter is conducted 
in a transverse direction toward the flow channels. The temperature is a maximum at 
any midpoint between adjacent channels, and the collector plate acts as a fin attached 
to the walls of the flow channel. The thermal performance of a fin can be expressed in 
terms of its efficiency. The fin efficiency r] f is defined as the ratio of the rate of heat 
flow through the real fin to the rate of heat flow through a fin of infinite thermal con¬ 
ductivity, that is, a fin at a uniform temperature. We shall now derive a relation to eval¬ 
uate this efficiency for a flat-plate solar collector. 

If U c is the overall unit conductance from the collector-plate surface to the ambi¬ 
ent air, the rate of heat loss from a given segment of the collector plate at x,y in Fig. 
3.14 is 


q(x,y) = U c [T c (x,y) - TJdxdy , (3.24) 

where T c = local collector-plate temperature (T c > T a ) 

T a = ambient air temperature 

U c = overall unit conductance between the plate and the ambient air 
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Figure 3.14. Sketch showing coordinates and dimensions for collector plate and fluid ducts. 
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U c includes the effects of radiation and free convection between the plates, the radia¬ 
tive and convective transfer between the top of the cover and the environment, and 
conduction through the insulation. Its quantitative evaluation has been previously con¬ 
sidered. 

If conduction in the x direction is negligible, a heat balance at a given distance x Q 
for a cross section of the flat-plate collector per unit length in the x direction can be 
written in the form 


a s I s dy - U C (T C - T a )dy + ( -kt^ 






y + dy,x o 


= 0. (3.25) 


If the plate thickness t is uniform and the thermal conductivity of the plate is inde¬ 
pendent of temperature, the last term in Eq. (3.25) is 


d Jc 

dy 


_ dT c 


+ 


ys o 


d%\ 
dy 2 ) 


dy. 


.Vv*o 


y+dys o dy 

and Eq. (3.25) can be cast into the form of a second-order differential equation: 

a± 


d% _ U c 
dy 2 kt ~ c 


T a + — 

a U, 


(3.26) 


The boundary conditions for the system described above at a fixed x (P are: 


1. At the center between any two ducts the heat flow is 0, or at y = 0, dT = 0. 

2. At the duct the plate temperature is T h (x 0 ), or at y = w = (/' - D)/2, V = T b (x 0 ), 
where T b (^ 0 ) is the fin-base temperature. 


If we let m 2 - UJkt and </) = 7 - (T a + a s IJU c ), Eq. (3.26) becomes 

d 2 $ 


— 


dy 


m 2 (f>. 


(3.27) 


subject to the boundary conditions 


defy 

—— = 0 at y = 0, and 
dy 

4> = T h (x 0 ) - (r a + aty = w. 


The general solution of Eq. (3.27) is 


0 = Cj sinh my + C 2 cosh my. 


(3.28) 
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The constants C, and C 0 can be determined by substituting the two boundary con¬ 
ditions and solving the two resulting equations for C, and C 2 . This gives 

lc - ( T a ± a JJ U c'> = cosh my 
T h (x 0 ) - (T a + ajjU c ) coshww' 

From the preceding equation the rate of heat transfer to the conduit from the por¬ 
tion of the plate between two conduits can be determined by evaluating the tempera¬ 
ture gradient at the base of the fin, or 

= — {aj s - U c [T b (x 0 ) - rjtanhmw}. (3.30) 

y=w m 

Since the conduit is connected to fins on both sides, the total rate of heat transfer is 


; dT c 

if. - -%r 


- u c [T b (x 0 ) - rj} 


tanh/mv 

mw 


(3.31) 


If the entire fin were at the temperature T b (jt 0 ), a situation corresponding physically to 
a plate of infinitely large thermal conductivity, the rate of heat transfer would be a 
maximum, # to[al max . As mentioned previously, the ratio of the rate of heat transfer with 
a real fin to the maximum rate obtainable is the fin efficiency 7] f . With this definition, 
Eq. (3.31) can be written in the form 


4W(*o) = 2w V f { a Js ~ U c lT b (x 0 ) ~T a \\, (3.32) 


where r] f = tanh mw/mw. 

The fin efficiency r) f is plotted as a function of the dimensionless parameter 
w(\J c /kt) ]/2 in Fig. 3.15. When the fin efficiency approaches unity, the maximum por¬ 
tion of the radiant energy impinging on the fin becomes available for heating the 
fluid. 

In addition to the heat transferred through the fin, the energy impinging on the 
portion of the plate above the flow passage is also useful. The rate of useful energy 
from this region available to heat the working fluid is 


= D{a; s - U c [T b (x 0 ) -T a ]}. (3.33) 

Thus, the useful energy per unit length in the flow direction becomes 

q u (x 0 ) = (D + 2w V )[aJf s - U c [T b (x 0 ) - TJ}. (3.34) 

The energy q u (x 0 ) must be transferred as heat to the working fluid. If the thermal 
resistance of the metal wall of the flow duct is negligibly small and there is no contact 
resistance between the duct and the plate, the rate of heat transfer to the fluid is 
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Figure 3.15. Fin efficiency for tube and sheet 
flat-plate solar collectors. 


q u (x 0 ) = Ph ci [T b (x 0 ) - T f (x 0 )l (3.35) 

where P is the perimeter of the flow duct, which is 2(D + d) for a rectangular duct. 
Contact resistance may become important in poorly manufactured collectors in which 
the flow duct is clamped or glued to the collector plate. Collectors manufactured by 
such methods are usually not satisfactory. 

3.2.6 Collector Efficiency Factor 

To obtain a relation for the useful energy delivered by a collector in terms of known 
physical parameters, the fluid temperature, and the ambient temperature, the collector 
temperature must be eliminated from Eqs. (3.34) and (3.35). Solving for T b (x 0 ) in Eq. 
(3.35) and substituting this relation in Eq. (3.34) gives 


q u (x 0 ) = l'F'{a s I s - U c [T f (x 0 ) - T a ]}, 


(3.36) 


where F' is called the collector efficiency factor [9], and V = (2w+D). F f is given by 


F' = 


_ l/Uc 

lil/(U e (D + 2wt ?/ » + \/(h c ,P)\ 


(3.37) 


For a flow duct of circular cross section of diameter D, F' can be written as 
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Table 3.6. Typical values for the parameters that determine the collector 
efficiency factor F f for a flat-plate collector in Eq. (3.37). 


U c 

2 glass covers 

4 W/m 2 • K 

0.685 Btu/hr • ft 2 • °F 

1 glass cover 

8 W/m 2 • K 

1.37 B Btu/hr • ft 2 • °F 

kt 

copper plate, 1 mm thick 

0.38 W/K 

0.72 Btu/hr • °F 

steel plate, 1 mm thick 

0.045 W/K 

0.0866 Btu/hr • °F 


water in laminar flow forced convection 

300 W/m 2 • K 

52 Btu/hr • ft 2 • °F 

water in turbulent flow forced convection 

1500 W/m 2 • K 

254 Btu/hr • ft 2 • °F 

air in turbulent forced convection 

100 W/m 2 • K 

17.6 Btu/hr • ft 2 • °F 


r= _ 1 /U e 

1'WKU c (D + 2wrj f )) + \/{h ci irD)\ 

Physically, the denominator in Eq. (3.37) is the thermal resistance between the fluid 
and the environment, whereas the numerator is the thermal resistance between the col¬ 
lector surface and the ambient air. The collector-plate efficiency factor F* depends on 
U c , h ci , and 7 ] f . It is only slightly dependent on temperature and can, for all practical 
purposes, be treated as a design parameter. Typical values for the factors determining 
the value of F' are given in Table 3.6. 

The collector efficiency factor increases with increasing plate thickness and plate 
thermal conductivity, but it decreases with increasing distance between flow channels. 
Also, increasing the heat-transfer coefficient between the walls of the flow channel 
and the working fluid increases F\ but an increase in the overall conductance U c will 
cause F' to decrease. 


3.2.7 Collector Heat-Removal Factor 

Equation (3.36) yields the rate of heat transfer to the working fluid at a given point x 
along the plate for specified collector and fluid temperatures. However, in a real col¬ 
lector the fluid temperature increases in the direction of flow as heat is transferred to it. 
An energy balance for a section of flow duct dx can be written in the form 

™ c P ( T f\x+dx - T f\x) = (3- 38 ) 

Substituting Eq. (3.36) for q u (x) and T^(x) 4- (dT f (x)/dx) dx for T^| x + dx in Eq. (3.38) 
gives the differential equation 
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mc p dT ^r = l ' F 'M - u c [T f (x) - rj}. 

Separating the variables gives, after some rearranging, 

_ dTfWi _ I'F'U 

T f (x) - T a - a s I s /U c mc p 


(3.39) 


(3.40) 


Equation (3.40) can be integrated and solved for the outlet temperature of the fluid 
T/.oiu f° r a duct l en g t h L, and for the fluid inlet temperature T /in if we assume that F 
and U c are constant, or 


Tlbul ~ T a- = ex J VJ'F'L 

T f,in~ T a~ a Js/ U c \ ™ C p 


(3.41) 


To compare the performance of a real collector with the thermodynamic optimum, it is 
convenient to define the heat-removal factor F R as the ratio between the actual rate of 
heat transfer to the working fluid and the rate of heat transfer at the minimum tempera¬ 
ture difference between the absorber and the environment. 

The thermodynamic limit corresponds to the condition of the working fluid re¬ 
maining at the inlet temperature throughout the collector. This can be approached 
when the fluid velocity is very high. From its definition F R can be expressed as 


Fr 


GCpiTfout T fin ) 

«Js - U c (T fJn - T a y 


(3.42) 


where G is the flow rate per unit surface area of collector m/A c . By regrouping the 
right-hand side of Eq. (3.42) and combining with Eq. (3.41), it can easily be verified 
that 


Fr 


Gc p j _ a .J s /U c (F f o Ut T a ) 
U c L a s I s /U c — (T f in — T a ) _ 




(3.43) 


Inspection of the above relation shows that F R increases with increasing flow rate and 
approaches as an upper limit F, the collector efficiency factor. Since the numerator of 
the right-hand side of Eq. (3.42) is q u , the rate of useful heat transfer can now be ex¬ 
pressed in terms of the fluid inlet temperature, or 


4u = A cFrWs - U c (Tf . n - T a )\. 

If a glazing above the absorber plate has transmittance r s then 


(3.44) 
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q u = A c F R [r s aJ c - U c (T fJn - T a )l (3.45) 

and instantaneous efficiency tj c is 

Vc = jj = Wa - U c (T Un - T a )/I c \ (3.46) 

Equation (3.46) is also known as the Hottel-Whillier-Bliss equation. This is a conve¬ 
nient form for design, because the fluid inlet temperature to the collector is usually 
known or can be specified. 


Example 3.1. Calculate the averaged hourly and daily efficiency of a water solar 
collector on January 15, in Boulder, CO. The collector is tilted at an angle of 60° and 
has an overall conductance of 8.0 W/m 2 • K on the upper surface. It is made of copper 
tubes, with a 1 -cm ID, 0.05 cm thick, which are connected by a 0.05-cm-thick plate at 
a center-to-center distance of 15 cm. The heat-transfer coefficient for the water in the 
tubes is 1500 W/m 2 • K, the cover transmittance is 0.9, and the solar absorptance of the 
copper surface is 0.9. The collector is 1 m wide and 2 m long, the water inlet tempera¬ 
ture is 330 K, and the water flow rate is 0.02 kg/sec. The horizontal insolation (total) l h 
and the environmental temperature are tabulated below. Assume the diffuse radiation 
accounts for 25 percent of the total insolation. 


Solution. The total radiation received by the collector is calculated from Eq. 
(2.42) and neglecting the ground reflected radiation: 

4 = Uc + h,c = 0-257, cos 2 (yj + (1 - 0.25)7,7?,. 


Time 

(hr) 

K 

(W/m 2 ) 

T 

amb 

(K) 

7-8 

12 

270 

8-9 

80 

280 

9-10 

192 

283 

10-11 

320 

286 

11-12 

460 

290 

12-13 

474 

290 

13-14 

395 

288 

14-15 

287 

288 

15-16 

141 

284 

16-17 

32 

280 


The tilt factor R b is obtained from its definition in Chapter 2 [see Eq. (2.58)]: 
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cos i _ sin(L - /3)sin 8 S H- cos(L — f3)co$8 s co sh s 
h sin a sin L sin 8 S 4- cos L cos6 v cos h s 

where L = 40° 8. = —21.1 on January 15 (from Fig. 2.8), and (3 = 60°. The hour angle 
/? s equals 15° for each hour away from noon. 

The fin efficiency is obtained from Eq. (3.32): 

_ tanh m( V - D)/2 
Vf ~ m{l' - D)/2 


where 


t/,V /2 

kt 


Vf 


390 X 5 X 10~ 4 
tanh6.4(0.15 - 0.01)/2 


1/2 


= 6.4, and 


= 0.938. 


6.4(0.15 - 0.01)/2 
The collector efficiency factor F' is, from Eq. (3.37), 

i/t4_ 


F' = 


l'[\/{U c {D + Iwr],)) + \/(h c ,TrD)\ 

1/8.0 


0.15[ 1/8.0(0.01 + 0.14 X 0.938) + 1/1500tt X 0.01] 
Then we obtain the heat-removal factor from Eq. (3.43): 


= 0.92. 


F„ = 


Gc, 

l/„ 


1 — exp - 


U C F' 
Gc„ 


Time 

(hr) 

h 

(W/m 2 ) 

R „ 

h, 

(W/m 2 ) 

4.c 

(W/m 2 ) 

/ 

(W/m 2 ) 

% 

(W) 

T 

amb 

(K) 

V c 

7-8 

12 

10.9 

1 

98 

99 

0 

270 

0 

8-9 

80 

3.22 

5 

193 

198 

0 

280 

0 

9-10 

192 

2.44 

12 

351 

363 

0 

283 

0 

10-11 

320 

2.18 

20 

523 

543 

148 

286 

0.137 

11-12 

460 

2.08 

29 

718 

747 

482 

290 

0.322 

12-13 

474 

2.08 

30 

739 

769 

512 

290 

0.333 

13-14 

395 

2.18 

25 

646 

671 

351 

288 

0.261 

14-15 

287 

2.49 

18 

525 

543 

175 

288 

0.162 

15-16 

141 

3.22 

9 

341 

350 

0 

284 

0 

16-17 

32 

10.9 

2 

261 

263 

0 

280 

0 
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Fr 


0.01 X 4184 

8.0 


1 — exp 


8.0 X 0,922 \ 
0.01 X 4184/ 


0.845. 


From Eq. (3.45), the useful heat delivery rate is 

q u = 2 X 0.845[/ c X 0.81 - 8.0(7},„ - T amb )). 

The efficiency of the collector is rj f = q u / AI c and the hourly averages are calculated 
in the table above. 

Thus, 2/ = 4546 W/m 2 and 2q u = 1668W. The daily average efficiency is ob¬ 
tained by summing the useful energy for those hours during which the collector deliv¬ 
ers heat and dividing by the total insolation between sunrise and sunset. This yields 

7 ] . = — = 0.183 or 18.3 percent 

hM' 2A r / c 2 X 4546 


3.2.8 Transient Effects 

The preceding analysis assumed that steady-state conditions exist during the operation of 
the collector. Under actual operating conditions the rate of insolation will vary and the 
ambient temperature and the external wind conditions may change. To determine the ef¬ 
fect of changes in these parameters on the performance of a collector it is necessary to 
make a transient analysis that takes the thermal capacity of the collector into account. 

As shown in [42], the effect of collector thermal capacitance is the sum of two 
contributions: the collector storage effect, resulting from the heat required to bring the 
collector up to its final operating temperature, and the transient effect, resulting from 
fluctuations in the meteorological conditions. Both effects result in a net loss of energy 
delivered compared with the predictions from the zero capacity analsysis. This loss is 
particularly important on a cold morning when all of the solar energy absorbed by the 
collector is used to heat the hardware and the working fluid, thus delaying the delivery 
of useful energy for some time after the sun has come up. 

Transient thermal analyses can be made with a high degree of precision [47], but 
the analytical predictions are no more accurate than the weather data and the overall 
collector conductance. For most engineering applications, a simpler approach is there¬ 
fore satisfactory [19]. For this approach, it will be assumed that the absorber plate, the 
ducts, the back insulation, and the working fluid are at the same temperature. If back 
losses are neglected, an energy balance on the collector plate and the working fluid for 
a single-glazed collector delivering no useful energy can be written in the form 

(SS),^jp = Va + A c U p [T s (t) - f p (t)l (3.47) 

where ( mc)p is the sum of the thermal capacities of the plate, the fluid, and the insola¬ 
tion, I s is the insolation on the absorber plate, and U p is the conductance between the 
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air-conditioning system. The magnitude of the reduction in heating system operation 
will be described in the next section. 


5.1.3 The Degree-day Method 

The preceding analysis of heat loss from buildings expresses the loss on a per unit temper¬ 
ature 1 ^ difference basis (except for unexposed floor slabs). In order to calculate the peak 
load and total annual load for a building, appropriate design temperatures mustbe define 
for each. The outdoor design temperature is usually defined statistically, such that the 
tual outdoor temperature will exceed the design temperature 97.5 or 99 percent ot the time 
over a long period. The design temperature difference (A7) is then the interior buildin 
temperature minus the outdoor design temperature. The design AT is used for rating no 
soli£ heating systems, but is not useful for selection of solar systems since solar systems 
rarely provide 100 percent of the energy demand of a building at peak conditions. 

A more useful index of heating energy demand is the total annual energy require¬ 
ment for a building. This quantity is somewhat more difficult to calculate than the pea 
load. It requires a knowledge of day-to-day variations in ambient temperature during 
the heating season and the corresponding building heat load for each day. Building 
heat loads vary with ambient temperatures as shown in Fig. 5.3. The environme 



Figure 5.3. Building load profile versus ambient temperature showing no-load 
temperature T, and desired interior temperature T r 
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temperature T nP above which no heat need be supplied to the building, is a few degrees 
below the required interior temperature T { because of internal heat-generation effects. 

The no-load temperature at which internal source generation q. just balances trans¬ 
mission and infiltration losses can be determined from the energy balance 

q, = UA{T { - T nl ), (5.5) 

where UA is the overall loss coefficient for the building (W/°C). Then 


T nl = T, - 


UA 


(5.6) 


The total annual heat load on the building, Q r can be expressed as 

Qr= f UA(T nl - T a ) + dt, (5.7) 

•'365 days 

in which all arguments of the integral are functions of time. The superscript + indicates 
that only positive values are considered. In practice, it is difficult to evaluate this inte¬ 
gral; therefore, three simplifying assumptions are made: 

1. UA is independent of time. 

2. T nl is independent of time. 

3. The integral can be expressed by the sum. 

Thus, 


_ 365 _ 

UA 2 {T ni - TX (5.8) 

n= 1 

where n is the day number, and the daily average temperature T a can be approximated 
by l A(T + T . ), in which T mov and T . are the daily maximum and minimum 

J v a, max a,m xn' 1 o.max a.min J 

temperatures, respectively. 

The quantity ( T nl — T a ) + is called the degree-day unit . For example, if the aver¬ 
age ambient temperature for a day is 5°C and the no-load temperature is 20°C, 15 de¬ 
gree C-days are said to exist for that day. However, if the ambient temperature is 20°C 
or higher, 0 degree-days exist, indicating 0 demand for heating that day. Degree-day 
totals for monthly (X month (T nl — T a ) + ) and annual periods can be used directly in Eq. 
(5.8) to calculate the monthly and annual heating energy requirements. 

In the past, a single value of temperature has been used throughout the United 
States as a universal degree-day base, 65.0°F or 18.3°C.* This practice is now out- 

*The degree-day base in SI units is defined as 19.0°C, not 18.3°C, which corresponds to 65.0°F. There¬ 
fore, precise conversion between the two systems is not possible by a simple multiplication by 5/9. 
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dated, since many homeowners and commercial building operators have lowered their 
thermostat settings in response to increased heating fuel costs, thereby lowering T nl . 
Likewise, warehouses and factories operate well below the 19°C level. Therefore, a 
more generalized database of degree-days to several bases (values of T nl ) has been cre¬ 
ated by the U.S. National Weather Service (NWS). 

Thom [32-34], in a series of papers, developed a statistically rigorous method of 
calculating degree-days to any base from values of ambient temperature and monthly 
standard deviations of ambient temperature. The details of this method are too lengthy 
to present here. However, the NWS has used this method to prepare tabulations of de¬ 
gree F-days for many U.S. locations. In addition, maps of degree F-days to the stan¬ 
dard base are available on a monthly basis [35]. 

A variable base degree-day method, which recognizes that T nl may vary not only 
with location but also from building to building, is more accepted now [3, 9, 12, 25]. 
Kreider and Rabl [25] have described this method in detail. Erbs et al. [12]_developed 
a model to estimate variable base degree-days that needs as input only T a for each 
month. Balcomb et al. [3] have listed monthly degree-days to bases of 50°F, 55°F, 
60°F, 65°F, and 70°F for 209 U.S. and Canadian cities. 


Example 5.2. A building located in Denver, CO, has a heat-loss coefficient UA of 
1000 kJ/hr • °C and internal heat sources of 4440 kJ/hr. If the interior temperature is 
20°C (68°F), what are the monthly and annual heating energy requirements? A gas fur¬ 
nace with 65 percent efficiency is used to heat the building. 


Solution. In order to determine the monthly degree-day totals, the no-load tem¬ 
perature (degree-day basis) must be evaluated from Eq. (5.6). 


Tnl = 20 - 


4440 

1000 


15.6° C(60° F) 


The monthly degree C-days for Denver are taken from the U.S. National Weather Ser¬ 
vice and given in Table 5.3. The energy demand is calculated as 

Energy Demand = UA X 24 — X Degree C — days. (5.9) 

day 

The monthly energy demand is given in Table 5.3. 

The annual energy demand of 62.9 GJ is delivered by a 65 percent efficient de¬ 
vice. Therefore, 

62.9 

Average Annual Purchased Energy = --^GJ = 96.8 GJ. 
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Table 5.3. Monthly and annual energy demands for Example 5.2. 


Month 


Degree C-days 

Energy demand 3 
(GJ) 

Jan. 


518 

12.4 

Feb. 


423 

10.2 

Mar. 


396 

9.5 

Apr. 


214 

5.2 

May 


68 

1.6 

June 


14 

0.3 

July 


0 

0 

Aug. 


0 

0 

Sep. 


26 

0.6 

Oct. 


148 

3.6 

Nov. 


343 

8.2 

Dec. 

Total 

472 

11.3 


2622 

62.9 


a Energy demand equals UA X degree C-days X 24 hrs/day. 


5.1.4 Service Hot-Water Load Calculation 

Service hot-water loads can be calculated precisely with the knowledge of only a few 
variables. The data required for calculation of hot-water demand are 

Water source temperature (T s ) 

Water delivery temperature ( T d ) 

Volumetric demand rate (Q) 

The energy requirement for service water heating q hw is given by 


RtJf) = P w Q(t)c pw [T d - T s (t)l (5.10) 

where p w is the water density and c is its specific heat. The demand rate, Q(t ), varies 
in general with time of day and time of year; likewise, the source temperature varies 
seasonally. Source temperature data are not compiled in a single reference; local water 
authorities are the source of such temperature data. 

Few generalized data exist with which to predict the demand rate Q. Table 5.4 in¬ 
dicates some typical usage rates for several common building types. Process water 
heating rates are peculiar to each process and can be ascertained by reference to pro¬ 
cess specifications. 

Example 5.3. Calculate the monthly energy required to heat water for a family of 
four in Nashville, TN. Monthly source temperatures for Nashville are shown in Table 
5.5, and the water delivery temperature is 60°C (140°F). 
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Table 5.4. Approximate service hot-water demand rates 



Demand per person 

Usage type 

liters/day 

gal/day 

Retail store 

2.8 

0.75 

Elementary school 

5.7 

1.5 

Multifamily residence 

76.0 

20.0 

Single-family residence 

76.0 

20.0 

Office building 

11.0 

3.0 


Solution. For a family of four, the demand rate Q may be found using a demand 
recommended from Table 5.4: 

Q — 4 X 76 liters/day = 0.30 m 3 /day 

The density of water can be taken as 1000 kg/m 3 and the specific heat as 4.18 kJ/kg • 
°C. 

Monthly demands are given by 

q m = (Q X days/month)(p H ,c /)lv )[7’ (/ - T s (t)] 

= (0.30 X days/month)( 1000 X 4.18)[60 - T s (t)]. 

The monthly energy demands calculated from the equation above with these data are 
tabulated in Table 5.5. 


Table 5.5. Water heating energy demands for Example 5.3. 


Month 

Days/month 

Demand 

(m-Vmonth) 

Source temperature 

C° 

Energy 

requirement 

(GJ/month) 

Jan. 

31 

9.3 

8 

2.0 

Feb. 

28 

8.4 

8 

1.8 

Mar. 

31 

9.3 

12 

1.9 

Apr. 

30 

9.0 

19 

1.5 

May 

31 

9.3 

17 

1.7 

June 

30 

9.0 

21 

1.5 

July 

31 

9.3 

22 

1.5 

Aug. 

31 

9.3 

24 

1.4 

Sep. 

30 

9.0 

24 

1.4 

Oct. 

31 

9.3 

22 

1.5 

Nov. 

30 

9.0 

14 

1.7 

Dec. 

31 

9.3 

12 

1.9 
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5.2 SOLAR WATER HEATING SYSTEMS 

Solar water-heating systems represent the most common application of solar energy at 
the present time. Small systems are used for domestic hot water applications while 
larger systems are used in industrial process heat applications. There are basically two 
types of water-heating systems: natural circulation or passive solar system (ther¬ 
mosyphon) and forced circulation or active solar system. Natural circulation solar wa¬ 
ter heaters are simple in design and low cost. Their application is usually limited to 
nonfreezing climates, although they may also be designed with heat exchangers for 
mild freezing climates. Forced circulation water heaters are used in freezing climates 
and for commercial and industrial process heat. 

5.2.1 Natural Circulation Systems 

The natural tendency of a less dense fluid to rise above a denser fluid can be used in a 
simple solar water heater to cause fluid motion through a collector [10]. The density 
difference is created within the solar collector where heat is added to the liquid. In the 
system shown in Fig. 5.4, as water gets heated in the collector, it rises to the tank, and 
the cooler water from the tank moves to the bottom of the collector, setting up a natural 
circulation loop. It is also called a thermosyphon loop. Since this water heater does not 
use a pump, it is a passive water heater. For the thermosyphon to work, the storage 
tank must be located higher than the collector. 

The flow pressure drop in the fluid loop (AP^^) must equal the bouyant force 
“pressure difference” (AP bouyant ) caused by the differing densities in the hot and cold 
legs of the fluid loop: 


Safety Valve Hot-water to Building 



Figure 5.4. Schematic diagram of thermosyphon loop used in a natural circulation, service water¬ 
heating system. The flow pressure drop in the fluid loop must equal the bouyant force “pressure” 

gp(x)dx - where p(x) is the local collector fluid density and p stor is the tank fluid density, 

assumed uniform. 
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AP, 


FLOW 


— AP BO uyant 



(5.11) 


where H is the height of the legs and L the height of the collector (see Fig. 5.4), p(x) is 
the local collector fluid density, p stor is the tank fluid density, and p oul is the collector 
outlet fluid density, the latter two densities assumed uniform. The flow pressure term 
APf Low’ > s related to the flow loop system head loss, which is in turn directly con¬ 
nected to friction and fitting losses and the loop flow rate: 


c1> pd(h,), 

Jloop 


(5.12) 


where h L = KV 2 , with K being the sum of the component loss velocity factors (see any 
fluid mechanics text) and V the flow velocity. 

Since the driving force in a thermosyphon system is only a small density differ¬ 
ence and not a pump, larger-than-normal plumbing fixtures must be used to reduce 
pipe friction losses [30]. In general, one pipe size larger than normal would be used 
with a pump system is satisfactory. Figure 5.5 shows some passive water heaters. 

Since the hot-water system loads vary little during a year, the angle of tilt is that 
equal to the latitude, that is, /3 = L. The temperature difference between the collector 
inlet water and the collector outlet water is usually 8—11°C during the middle of a 
sunny day [10]. After sunset, a thermosyphon system can reverse its flow direction and 
lose heat to the environment during the night. To avoid reverse flow, the top header of 
the absorber should be at least 30 cm below the cold leg fitting on the storage tank, as 
shown, otherwise a check valve would be needed. 

To provide heat during long cloudy periods, an electrical immersion heater can be 
used as a backup for the solar system. The immersion heater is located near the top of the 
tank to enhance stratification and so that the heated fluid is at the required delivery tem¬ 
perature. Tank stratification is desirable in a thermosyphon to maintain flow rates as high 
as possible. Insulation must be applied over the entire tank surface to control heat loss. 

Several features inherent in the thermosyphon design limit its utility. If it is to be 
operated in a freezing climate, a nonfreezing fluid must be used, which in turn requires 
a heat exchanger between collector and potable water storage. (If potable water is not 
required, the collector can be drained during cold periods instead.) Heat exchangers of 
either the shell-and-tube type or the immersion-coil type require higher flow rates for 
efficient operation than a thermosyphon can provide. Therefore, the thermosyphon is 
usually limited to nonfreezing climates. For mild freeze climates, a heat exchanger coil 
welded to the outer surface of the tank and filled with an antifreeze may work well. 

Example 5.4. Determine the “pressure difference” available for a thermosyphon 
system with 1 meter high collector and 2 meter high legs. The water temperature input 
to the collector is 25°C and the collector output temperature is 35°C. If the overall sys¬ 
tem loss velocity factor (K) is 15.6, estimate the system flow velocity. 
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Figure 5.5. Passive solar water heaters; (a) compact model using combined collector and storage, 
(b) section view of the compact model, and (c) tank and collector assembly. 


Solution. Equation (5.11) is used to calculate the pressure difference, with the 
water densities being found from the steam tables (see Appendix 3 and Tables 3.8 and 
3.9). 


p stor (25°C) = 997.009 kg/m 3 
P OU t(35°C) = 994.036 kg/m 3 

Pcoii.ave.(30°O = 996.016 kg/m 3 (note; average collector temperature used in ‘inte¬ 
gral’) and with H = 2 and L = 1 m: 

APbouvant = (997.009)9.81(2) - [(996.016)9.81(1) + (994.036)9.81(1)] 

= 38.9 N/m 2 (Pa). 

The system flow velocity is estimated from the system K given, the pressure difference 
calculated above, taking the average density of the water around the loop (at 30°C), 
and substituting into Eq. (5.12); 
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^^BOUYANT (Ploop.ave.X^i/loop (Ploop.ave.) 

V 2 = 38.9/(996.016)(15.6) 

V = 0.05 m/s 

5.2.2 Forced-Circulation Systems 

If a thermosyphon system cannot be used for climatic, structural, or architectural rea¬ 
sons, a forced-circulation system is required. 

Figure 5.6 shows three configurations of forced circulation systems: (1) open 
loop, (2) closed loop, and (3) closed loop with drainback. In an open loop system (Fig. 
5.6a) the solar loop is at atmospheric pressure, therefore, the collectors are empty 
when they are not providing useful heat. A disadvantage of this system is the high 
pumping power required to pump the water to the collectors every time the collectors 
become hot. This disadvantage is overcome in the pressurized closed loop system (Fig. 
5.6b) since the pump has to overcome only the resistance of the pipes. In this system, 
the solar loop remains filled with water under pressure. 

In order to accommodate the thermal expansion of water from heating, a small 
(about 2 gallon capacity) expansion tank and a pressure relief valve are provided in the 
solar loop. Because water always stays in the collectors of this system, antifreeze 
(propylene glycol or ethylene glycol) is required for locations where freezing condi¬ 
tions can occur. During stagnation conditions (in summer), the temperature in the col- 



Figure 5.6. Typical configurations of solar water-heating systems: (a) open loop system. 


{continues) 
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(b) 



Figure 5.6. ( Continued) (b) closed loop system; and (c) closed loop drainback system. 
(Adapted from Goswami, D.Y. Alternative Energy in Agriculture, Vol. 1. Boca Raton, FL: 
CRC Press, (1986). 




































SOLAR-HEATING SYSTEMS 


227 


lector can become very high, causing the pressure in the loop to increase. This can 
cause leaks in the loop unless some fluid is allowed to escape through a pressure- 
release valve. Whether the result of leaks or of draining, air enters the loop causing the 
pumps to run dry. This disadvantage can be overcome in a closed loop drainback sys¬ 
tem which is not pressurized (Fig. 5.6c). In this system, when the pump shuts off, the 
water in the collectors drains back into a small holding tank while the air in the holding 
tank goes up to fill the collectors. The holding tank can be located where freezing does 
not occur, but still at a high level to reduce pumping power. In all three configurations, 
a differential controller measures the temperature differential between the solar collec¬ 
tor and the storage, and turns the circulation pump on when the differential is more 
than a set limit (usually 5°C) and turns it off when the differential goes below a set 
limit (usually 2°C). Alternatively, a photovoltaic (PV) panel and a DC pump may be 
used. The PV panel will turn on the pump only when solar radiation is above a mini¬ 
mum level. Therefore, the differential controller and the temperature sensors may be 
eliminated. 


5.2.3 Industrial Process Heat Systems 

For temperatures of up to about 100°C, required for many industrial process heat ap¬ 
plications, forced emulation water-heating systems described above can be used. The 
systems, however, will require a large collector area, storage and pumps, etc. For 
higher temperatures, evacuated tube collectors or concentrating collectors must be 
used. Industrial process heat systems are described in more detail in Chapter 8. 


5.3 LIQUID-BASED SOLAR-HEATING SYSTEMS FOR BUILDINGS 

The earliest active solar space-heating systems were constructed from enlarged water¬ 
heating components. Experiments beginning in 1938 at the Massachusetts Institute of 
Technology (MIT) showed that solar heating with liquid working fluids could be done 
without any major technical problems. The early MIT work formed the basis of many 
of the design techniques used today. Other experiments after World War II provided 
additional fundamental information on collector designs and storage operation for liq¬ 
uid-based heating systems. Fig. 5.7 shows a modern solar-heated, multitamily building 
in Colorado. 

Solar space-heating systems can be classified as active or passive depending on 
the method utilized for heat transfer. A system that uses pumps and/or blowers for fluid 
flow in order to transfer heat is called an active system. On the other hand, a system 
that utilizes natural phenomena for heat transfer is called a passive system. Examples 
of passive solar space-heating systems include direct gain, attached greenhouse, and 
storage wall (also called Trombe wall). Passive solar heating systems are described in 
Chapter 7. In this section, configurations, design methods, and control strategies for 
active solar-heating systems are described. 
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Figure 5.7. Multi-unit residence located in Boulder, CO, heated by an active solar¬ 
heating system. Courtesy of Joint Venture, Inc., Boulder, CO. 


5.3.1 Physical Configurations of Active Solar Heating Systems 

Figure 5.8 is a schematic diagram of a typical space-heating system. The system con¬ 
sists of three fluid loops-collector, storage, and load. In addition, most space-heating 
systems are integrated with a domestic water-heating system to improve the year long 
solar load factor. 

Since space heating is a relatively low-temperature use of solar energy, a thermo¬ 
dynamic match of collector to task indicates that an efficient flat-plate collector or 
low-concentration solar collector is the thermal device of choice. 

The collector fluid loop contains fluid manifolds, the collectors, the collector 
pump and heat exchanger, an expansion tank, and other subsidiary components. A col¬ 
lector heat-exchanger and antifreeze in the collector loop are normally used in all solar 
space-heating systems, since the existence of a significant heating demand implies the 
existence of some subfreezing weather. 

The storage loop contains the storage tank and pump as well as the tube side of the 
collector heat exchanger. To capitalize on whatever stratification may exist in the stor¬ 
age tank, fluid entering the collector heat exchanger is generally removed from the 
bottom of storage. This strategy ensures that the lowest temperature fluid available in 
the collector loop is introduced at the collector inlet for high efficiency. 

The energy delivery-to-load loop contains the load device, baseboard heaters or 
fin-and-tube coils, and the backup system with a flow control (mode selector) valve. 

5.3.2 Solar Collector Orientation 

The best solar collector orientation is such that the average solar incidence angle is 
smallest during the heating season. For tracking collectors this objective is automati- 




Air bleed Relief valve Mode selector 
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Figure 5.8. Typical solar-thermal system for space heating and hot-water heating showing fluid transport loops and pumps. 






















































Table A2.6a Solar position and insolation values for 24 degrees north latitude 0 



Solar 

Solar 








Date 

time 

position 


BTUH/sq. ft. total insolation on surface* 










South facing surface 



AM 

PM 

Alt 

Azm 




angle with horiz. 







Normal 0 

Horiz. 

14 

24 

34 

44 

90 

Jan 21 

7 

5 

4.8 

65.6 

71 

10 

17 

21 

25 

28 

31 


8 

4 

16.9 

58.3 

239 

83 

110 

126 

137 

145 

127 


9 

3 

27.9 

48.8 

288 

151 

188 

207 

221 

228 

176 


10 

2 

37.2 

36.1 

308 

204 

246 

268 

282 

287 

207 


11 

1 

43.6 

19.6 

317 

237 

283 

306 

319 

324 

226 


12 

46.0 

0.0 

320 

249 

296 

319 

332 

336 

232 



Surface daily totals 

2766 

1622 

1984 

2174 

2300 

2360 

1766 

Feb 21 

7 

5 

9.3 

74.6 

158 

35 

44 

49 

53 

56 

46 


8 

4 

22.3 

67.2 

263 

116 

135 

145 

150 

151 

102 


9 

3 

34.4 

57.6 

298 

187 

213 

225 

230 

228 

141 


10 

2 

45.1 

44.2 

314 

241 

273 

286 

291 

287 

168 


11 

1 

53.0 

25.0 

321 

276 

310 

324 

328 

323 

185 


12 

56.0 

0.0 

324 

288 

323 

337 

341 

335 

191 



Surface daily totals 

3036 

1998 

2276 

2396 

2436 

2424 

1476 

Mar 21 

7 

5 

13.7 

83.3 

194 

60 

63 

64 

62 

59 

27 


8 

4 

27.2 

76.8 

267 

141 

150 

152 

149 

142 

64 


9 

3 

40.2 

67.9 

295 

212 

226 

229 

225 

214 

95 


10 

2 

52.3 

54.8 

309 

266 

285 

288 

283 

270 

120 


11 

1 

61.9 

33.4 

315 

300 

322 

326 

320 

305 

135 


12 

66.0 

0.0 

317 

312 

334 

339 

333 

317 

140 



Surface daily totals 

3078 

2270 

2428 

2456 

2412 

2298 

1022 

Apr 21 

6 

6 

4.7 

100.6 

40 

7 

5 

4 

4 

3 

2 


7 

5 

18.3 

94.9 

203 

83 

77 

70 

62 

51 

10 


8 

4 

32.0 

89.0 

256 

160 

157 

149 

137 

122 

16 


9 

3 

45.6 

81.9 

280 

227 

227 

220 

206 

186 

46 


10 

2 

59.0 

71.8 

292 

278 

282 

275 

259 

237 

61 


11 

1 

71.1 

51.6 

298 

310 

316 

309 

293 

269 

74 


12 

77.6 

0.0 

299 

321 

328 

321 

305 

280 

79 



Surface daily totals 

3036 

2454 

2458 

2374 

2228 

2016 

488 

May 21 

6 

6 

8.0 

108.4 

86 

22 

15 

10 

9 

9 

5 


7 

5 

21.2 

103.2 

203 

98 

85 

73 

59 

44 

12 


8 

4 

34.6 

98.5 

248 

171 

159 

145 

127 

106 

15 


9 

3 

48.3 

93.6 

269 

233 

224 

210 

190 

165 

16 


10 

2 

62.0 

87.7 

280 

281 

275 

261 

239 

211 

22 


11 

1 

75.5 

76.9 

286 

311 

307 

293 

270 

240 

34 


12 

86.0 

0.0 

288 

322 

317 

304 

281 

250 

37 



Surface daily totals 

3032 

2556 

2447 

2286 

2072 

1800 

246 

Jun 21 

6 

6 

9.3 

111.6 

97 

29 

20 

12 

12 

11 

7 


7 

5 

22.3 

106.8 

201 

103 

87 

73 

58 

41 

13 


8 

4 

35.5 

102.6 

242 

173 

158 

142 

122 

99 

16 


9 

3 

49.0 

98.7 

263 

234 

221 

204 

182 

155 

18 


10 

2 

62.6 

95.0 

274 

280 

269 

253 

229 

199 

18 


11 

1 

76.3 

90.8 

279 

309 

300 

283 

259 

227 

19 


12 

89.4 

0.0 

281 

319 

310 

294 

269 

236 

22 


Surface daily totals 

2994 

2574 

2422 

2230 

1992 

1700 

204 


590 


























































Table A2.6a Solar position and insolation values for 24 degrees north 
latitude" (Continued) 



Solar 

Solar 








Date 

time 

position 


BTUH/sq. ft. total insolation on surface^ 










South facing surface 



AM 

PM 

Alt 

Azm 




angle with horiz. 







Normal* 

Horiz. 

14 

24 

34 

44 

90 

Jul 21 

6 

6 

8.2 

109.0 

81 

23 

16 

11 

10 

9 

6 


7 

5 

21.4 

103.8 

195 

98 

85 

73 

59 

44 

13 


8 

4 

34.8 

99.2 

239 

169 

157 

143 

125 

104 

16 


9 

3 

48.4 

94.5 

261 

231 

221 

207 

187 

161 

18 


10 

2 

62.1 

89.0 

272 

278 

270 

256 

235 

206 

21 


11 

1 

75.7 

79.2 

278 

307 

302 

287 

265 

235 

32 


12 

86.6 

0.0 

280 

317 

312 

298 

275 

245 

36 



Surface daily totals 

2932 

2526 

2412 

2250 

2036 

1766 

246 

Aug 21 

6 

6 

5.0 

101.3 

35 

7 

5 

4 

4 

4 

2 

7 

5 

18.5 

95.6 

186 

82 

76 

69 

60 

50 

11 


8 

4 

32.2 

89.7 

241 

158 

154 

146 

134 

118 

16 


9 

3 

45.9 

82.9 

265 

223 

222 

214 

200 

181 

39 


10 

2 

59.3 

73.0 

278 

273 

275 

268 

252 

230 

58 


11 

1 

71.6 

53.2 

284 

304 

309 

301 

285 

261 

71 


12 

78.3 

0.0 

286 

315 

320 

313 

296 

272 

75 



Surface daily totals 

2864 

2408 

2402 

2316 

2168 

1958 

470 

Sep 21 

7 

5 

13.7 

83.8 

173 

57 

60 

60 

59 

56 

26 


8 

4 

27.2 

76.8 

248 

136 

144 

146 

143 

136 

62 


9 

3 

40.2 

67.9 

278 

205 

218 

221 

217 

206 

93 


10 

2 

52.3 

54.8 

292 

258 

275 

278 

273 

261 

116 


11 

1 

61.9 

33.4 

299 

291 

311 

315 

309 

295 

131 


12 

66.0 

0.0 

301 

302 

323 

327 

321 

306 

136 



Surface daily totals 

2878 

2194 

2342 

2366 

2322 

2212 

992 

Oct 21 

7 

5 

9.1 

74.1 

138 

32 

40 

45 

48 

50 

42 


8 

4 

22.0 

66.7 

247 

111 

129 

139 

144 

145 

99 


9 

3 

34.1 

57.1 

284 

180 

206 

217 

223 

221 

138 


10 

2 

44.7 

43.8 

301 

234 

265 

277 

282 

279 

165 


11 

1 

52.5 

24.7 

309 

268 

301 

315 

319 

314 

182 


12 

55.5 

0.0 

311 

279 

314 

328 

332 

327 

188 



Surface daily totals 

2868 

1928 

2198 

2314 

2364 

2346 

1442 

Nov 21 

7 

5 

4.9 

65.8 

67 

10 

16 

20 

24 

27 

29 


8 

4 

17.0 

58.4 

232 

82 

108 

123 

135 

142 

124 


9 

3 

28.0 

48.9 

282 

150 

186 

205 

217 

224 

172 


10 

2 

37.3 

36.3 

303 

203 

244 

265 

278 

283 

204 


11 

1 

43.8 

19.7 

312 

236 

280 

302 

316 

320 

222 


12 

46.2 

0.0 

315 

247 

293 

315 

328 

332 

228 



Surface daily totals 

2706 

1610 

1962 

2146 

2268 

2324 

1730 

Dec 21 

7 

5 

3.2 

62.6 

30 

3 

7 

9 

11 

12 

14 


8 

4 

14.9 

55.3 

225 

71 

99 

116 

129 

139 

130 


9 

3 

25.5 

46.0 

281 

137 

176 

198 

214 

223 

184 


10 

2 

34.3 

33.7 

304 

189 

234 

258 

275 

283 

217 


11 

1 

40.4 

18.2 

314 

221 

270 

295 

312 

320 

236 


12 

42.6 

0.0 

317 

232 

282 

308 

325 

332 

243 



Surface daily totals 

2624 

1474 

1852 

2058 

2204 

2286 

1808 


"From Kreider, J. F., and F. Kreith, “Solar Heating and Cooling,” revised 1st ed.. Hemisphere Publ. 
Corp., 1977. 

b \ Btu/hr • ft 2 = 3.152 W/m 2 . Ground reflection not included on normal or horizontal surfaces. 
^Normal insolation does not include diffuse component. 
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